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Abstract 

The contribution to the low frequency internal friction and the thermal conductivity due to 
optically vibrating edge dislocation dipoles is calculated within the modified Granato-Liicke string 
model. The results are compared with the recent experiments on plastically deformed samples 
of Al, Ta and Nb at low temperatures. It is shown that the presence of a reasonable density of 
optically vibrating dislocation dipoles provides a good fit to the thermal conductivity in super- 
conducting samples. At the same time, the internal friction experiments can not be described 
within the standard fluttering string mechanism. We found that the problem can be solved by 
assuming random forces acting on the dislocation dipoles. This gives an additional contribution 
to the internal friction which describes well the experimental data at low temperatures while their 
contribution to the thermal conductivity is found to be negligible. 

PACS numbers: 61.72.Lk, 62.40.+i, 61.72.Hh 
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I. INTRODUCTION 



Recently , 1 ^ 3 ^ 5 the low-temperature internal friction, thermal conductivity, specific heat, 
and heat release of plastically deformed, high-purity superconducting crystalline samples 
of Al, Nb, and Ta have been experimentally studied and compared with measurements on 
amorphous SiC>2 specimen.— In particular, it was established that plastic deformation 
has a pronounced effect on the internal friction and the thermal conductivity. Namely, the 
value of the internal friction can be increased by two orders of magnitude over that observed 
on annealed samples and it becomes comparable to that of amorphous SiC>2. Likewise, the 
thermal conductivity was found to have the similar value as that of amorphous SiCV-^ 
At the same time, neither long-time heat release nor additional contribution to the heat 
capacity expected for amorphous systems was observed. This finding indicates that the 
phonon scattering by dislocations appearing in crystals under plastic deformations could be 
of importance. 

As is well known, there are two principal mechanisms of phonon scattering by dislocations: 
the static strain-field scattering and the reradiation scattering.^ The first one is due to the 
anharmonicity of the dislocation strain field while the second one results from a possibility 
of the sound wave to induce a dislocation vibration. In this case, the incident energy will 
be dissipated as long as the dislocation radiates elastic waves. This is so-called fluttering 
mechanism which is dynamical in its origin. It is a difficult problem to experimentally 
clarify what kind of scattering actually dominates and especially to obtain a quantitative 
agreement between the observed effect and theoretical calculations based on either static 
or dynamic models. Moreover, it was shown that both mechanisms of phonon scattering 
by an array of single dislocations are failed in the description of various experiments.— i 3 i 4 i 8 i 9 
For example, the estimated contribution to the thermal conductivity due to the resonant 
interaction was found to have a reasonable agreement with the experimental data only at 
very short lengths of dislocation lines (otherwise the resonant frequency becomes too low). 
Notice that a similar problem has already emerged in experiments with plastically deformed 
LiF Namely, the measurements of the thermal conductivity and the ballistic phonon 
propagation in deformed LiF at low temperatures show that the obtained phonon scattering 
is too strong to be explained by static mechanisms of phonon-dislocation interaction, but 
is in rough agreement with calculations based on a resonant or dynamic interaction with 
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dislocations which can flutter in the stress field of passing phonons. The crucial role in 
explanation of the experimental data plays an assumption given in Reffiol that a reasonable 
density of optically vibrating dislocation dipoles is present in LiF. In this case, the resonant 
frequency becomes markedly higher. It should be mentioned that dislocation dipoles with 
a rather high density were actually observed in LiF— as well as in some other materials. 
For example, estimates from experiments of the ratio of the dislocation dipole density to a 
density of single dislocations range from 1000 (in strain-hardening measurements^) to 100 
(by deformation- induced bulk-density changes^) and less (by electron microscopyi 3 -^) . 

The purpose of this paper is to show that the concept of dynamical scattering of phonons 
by dislocation dipoles 1 ^ previously used for the explanation of the experiments in LiF&£ 
could also be suitable for Al, Nb, and Ta. As was observed in Refsjsfsf^ , these metals show 
a complex dislocation structure under plastic deformation. We assume that, by analogy with 
LiF, dislocation dipoles are also present in plastically deformed Al, Nb, and Ta. It should 
be noted that the dipole is a stable dislocation structure arising under plastic deformation 
in metals (see, e.g., Ref il5T) . Recently , the substructure of the dislocation ensemble in 
Al and Cu was experimentally studied under cyclic deformation. In particular, it was found 
in Ref r] that for Al deformed at 77 K the microstructure is composed exclusively of vein 
dipole bundles and channels without persistent slip bands. The dislocation density in the 
veins is about 3.2 x 10 14 m -2 , and the density in the channels is about 2.4 x 10 13 m -2 with 
the dipole separation in a range from 3 nm to 30 nm. It is important to note that all of 
the channel defects were found to be edge dipoles. In addition, the computer modelling 
performed in Refill shows that there is a critical size of the dipole separation (1.6 nm for 
Al and 0.42 nm for Cu) above which edge-dislocation dipoles in Al and Cu become stable 
with respect to the athermal annihilation. Notice that we consider in our paper dislocation 
dipoles in Al with a separation close to 20 nm, that is much larger the critical size found in 
Reffl 

Thus, as for LiF the gain in resonant frequency can be achieved due to optically vibrating 
dipoles having a reasonable dipole separation and dislocation length. This allows us to obtain 
a quantitative agreement with the experimental data for all three metals. Evidently, the 
presence of the fluttering dislocation dipoles will affect not only the thermal conductivity 
but also other physical characteristics, in particular, the internal friction. We study the 
contribution to the internal friction due to vibrating dislocation dipoles in the framework of 
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the modified Granato-Liicke string theory. As is known, the Granato-Liicke vibrating string 
model is based on an analogy between the vibration of a pinned dislocation line segment 
and the forced damped vibration of a string. 19 

Our study shows that a good agreement with the internal friction experiments in Al, Nb, 
and T a 1 ^ 5 can be obtained only in the presence of some random component of the stresses. 
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This idea has been recently proposed and studied in Refs|20| and |21[ In particular, it was 
found that the action of random forces on the dislocation gives a substantial contribution 
to the decrement. As was observed experimentally,—^ the decrement has an unexpectedly 
large value and practically does not depend on the amplitude and the frequency of external 
stresses. This means that the noise parameters should be governed by the characteristics 
of an external signal because this is the only way to provide the constant contribution to 
the decrement. Therefore, we suggest that the external force introduces also a specific time 
scale for correlations of stress fluctuations in dislocation ensemble which always has its own 
stochastic dynamics.— This time scale determines the time of an additional relaxation of 
internal stresses due to random stresses. It seems reasonable to assume that this time is 
related to the frequency of the external force. As the simplest approximation, we use the 
correlation function exponentially decreasing with time. What is important, the introduction 
of a time scale means that dislocations experience the influence of the colored noise. Finally, 
we adapt the theory proposed in RefsQ and to the case of dislocation dipoles and 
introduce the additional noise component with the above-mentioned correlation function in 
the equations of motion. 



II. THE MODEL 



Let us consider a dislocation dipole with dislocation lines directed along the z-axis and 
lying in parallel glide planes at distance I. The Burgers vectors are opposite while the 
orientation of dislocation lines is the same. The xz plane is chosen to be a slip plane for a 
positive dislocation. In this case the plane y = I is the slip plane for a negative dislocation. 
Within the string approximation, the dislocation dipole can be modeled by means of two 
damped interacting vibrating strings. Supposing that all the parameters of the string are 
the same for both dislocations, one can formulate the general equations of damped glide 
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motion for the edge dislocation dipole in the form 



de(z,t) d 2 e(z,t) t+ 
m <T7o r rf — 1q — — = t , 



_ &rj>(z,t) D m*,t) r gV(M) pe , t _ m 



<9t 2 ^ " «9,2 2 

>(z,t) m*,t) d 2 j;(z 

dt 2 dt ° dz 2 

where e(z,t) and 4>(z,t) are the displacements of the positive and negative dislocations, 
m the effective mass, T the line tension, B the damping parameter, and F ext± the total 
external force which acts on the positive and negative dislocations, correspondingly, in their 
glide planes. The parameters of the model are written as 1 - 

m=?f[l + (c t /c l )*]\ng, T = f \ng, B = ^[1 + (c t / Cl )% (2) 
An 47T (1 — v) 8 

where p is the density, p the shear modulus, ID the thermal phonon frequency, v the Poisson 
constant, b the length of Burgers vector, q and q the transverse and longitudinal sound 
velocities. At low frequencies, g can be estimated as g ~ 1/b^/X where A is the dislocation 
dipole density. The total external force, resolved in the glide plane, includes three terms 



f ± + F ± + b ± V (t). (3) 



Here are the interaction forces between the dislocations in the dipole, F^ 1 are the forces 
due to external stress field cr^, and rj(t) describes a stationary random component. The 
explicit form of r](t) is not specified at this stage. Hereafter, the sign plus (minus) corresponds 
to a positive (negative) dislocation, and the summation over repeated indices is assumed. 
In general, the interaction between dislocations is determined by the Peach-Koehler force 

ff = e jak r a bfaf k , (4) 

where Ej a k is the totally antisymmetric tensor, af k are the corresponding stresses of dislo- 
cations in the dipole, and f is the unit tangent vector to the defect line. The displacement 
fields can be written as 

u n (r,t) = J c ijkl G jnii 5e p k \dV' , (5) 

with Cyki being the elastic modulus, Gj n the Green's tensor, Gj n ^ = dGj n /dxi, and 5e p k \ the 
variation of the plastic part of the strain tensor. For sliding dislocation one has^ 

Sei = lihidxr^ + 6,[*ir|*W3. ( 6 ) 
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where 5x describes the displacement of the dislocation line from a straight configuration 
and <5(£) is the two-dimensional delta-function. In our geometry, b + = —b~ = (6,0,0), 
t = (0,0, —1), £ + = (x',y',0), and £~ = (x',y' — 1,0). The displacements are chosen to be 
5x+ = (e(z),0,0) and 5x~ = (ip(z), 0,0), so that [5x+t\ = (0, e, 0) and [5x~r\ = (0,^,0). 
After substitution into Eq.fjSJ) one obtains 



u t. (r) = J Cijk2G jn ,i{x, y,z- z')b%e(z')dz', 
u n(r) = / c ijk2 G jnti {x, y-l,z- z')b^i/)(z')dz' . 



(7) 



For isotropic case, 



G km (r) 



h 1 

r km 2(1 - v) 



I '.km 



and Cijki = XSijSki + ^(SikSji + 6u5jk) with A and \i being the Lame constants. Supposing 
that e(z) = e(k) exp(ikz) and ip(z) = i/j(k) exp(ikz), after straightforward calculations one 
obtains from Eq.(J3J 



/-(f) = -2e{k)e lkz ju 2 b 2 k 2 



2 A [ 1 - K 2 (ka) - AAkx y K^ka) + (2 A + B)K (ka) 



a. 



a" 



(8) 



/+(f) = -2^(k)e lkz fi 2 b 2 k 2 
+ (2A + B)K (kc) 



(9) 



where A = —5/4(1 — v), B = 1/Att/j,, a 2 = x 2 + y 2 , c 2 = x 2 + (y — I) 2 , and are the 
McDonalds functions. Notice that we are interested only in a component f\ since, at chosen 
geometry, all the remaining components become equal to zero in slip planes of dislocations. 
In the long wavelength approximation (hi <C 1) one gets 



f-{0,l,z) = De{z), f + (0,0,z) = D^(z), 



where 



D 



fib 2 



(10) 



(11) 



2tt(1 -v)l 2 ' 

The second term in Eq. © is also determined by the Peach-Koehler force which is written 
in the known form F^ 1 = s ra kT a bfaik with being an external stress field. In our geometry, 



F + = b + a(0,0,z), F~ — b~a(0,l,z). 



(12) 



Let us consider first the case of free vibrations of the dipole when B = F ± = rj(t) = in 
Eq.([Tj). There exist both acoustical and optical modes with e = ip and e = —ip, respectively. 
Substituting Eq.flTUJ) into Eq.(JTJ) one obtains 

d 2 e{z,t) d 2 e{z,t) 
m —fa2 T ° dz 2 ±De ^ = - ( 13 ) 

Here the sign plus corresponds to the optical mode while the minus to the acoustical one. 
The spectrum of normal oscillations is found to be 



n ^=r" ± ^' (14) 

where u n = {im / L) J 'Tq / m is the resonance frequency of a single dislocation and L is the 
dislocation length. For optical mode, this spectrum was presented in Ref. As is seen, the 
resonance frequency becomes higher for the optical mode and the difference increases with 
decreasing dipole separation. On the contrary, for acoustical mode the resonance frequency 
lies even lower than that for a single dislocation. 

The measured thermal conductivity can be fitted by assuming dislocation dipoles having 
a higher resonant frequency and being more numerous than isolated dislocations. 10 For this 
reason, we consider below only the case of optically vibrating dislocation dipoles. Supposing 
e = —if), cr(0, 0, z) = cr(0, 1, z), and taking into account Eqs.(J3J), (fTUj) and (fT2j) one can reduce 
Eq.© to 

d 2 e(z,t) m d 2 e(z, t) de(zA) , „ . ^ , . , .. . . 

m ^ } = T ^ } - + ba(0, 0, z) - De{z, t) + bq{t). (15) 

Let us consider a periodic external stress wave in the form 

a(0, 0, z) = roe™ = £ a n sm(k n z)e™ } (16) 

n 

where o"o is the shear stress component resolved in the glide plane, k n = nn/L, and a n = 
4ao/irn is the Fourier coefficient. Then the general solution to Eq. (fTH|) is written as 



e(z,t) = e n(t) sm(nnz/L). (17) 

n=l 

Besides, e n (t) can be presented as the sum of periodic and random components 

e n (t) = e n (t) p + e n (t) r . (18) 
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Substitution of Eqs.lfTTj) and (fTHj) into Eq.(JT3J) gives the following system of equations: 

d 2 e n (t) r B de n (t) r , 4^(t) _ 

at 1 m ot ixnm 

at 2 m at imm 
The solution to the second equation is written as 

e n (t) P = C n e™\ (20) 

with 

^ bcr n 1 



\ m n 

The solution to the first equation in Eq. (JT§j) is written in the form 



t 



"n 4m 2 

As is seen, the displacement e n (z,t) depends on the random function rj(t). As a result, the 
dipole dynamics is governed by a random function having its own probabilistic character- 
istics. This is an essential difference from the standard string model where the string is 
subjected only to the influence of the periodic force. 

III. INTERNAL FRICTION 

Let us calculate the total decrement A t due to the motion of optically vibrating edge 
dislocation dipoles by using the formulated model. Since e n {t) is a random function, one has 
to perform an additional averaging over the possible realizations of random stresses r](t). 
The mean energy loss to friction due to optically vibrating edge dislocation dipole per unit 
time reads 

L T . . 

<P>=2| T liml< J(F+ + ^(t))^±dt>dz, (22) 
o o 

where <> means averaging over random force realizations ensemble. Generally, 

__ < AW > _ N < AW d > _ N < P > 2tt 
t ~2<W>~ 2<W> ~2<W>~' ( ' 
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where < W > is the averaged total vibration energy per unit volume, N is the total number of 
dislocation dipoles per unit volume, < AWd > is the contribution from the single dislocation 
dipole. After substitution of Eq.(|18|) into Eq.(|22j) the decrement is written as 

A t = A r + A p , (24) 

where 

^ = & g J < ^ T - S ^ > e - BS/2m - £r^) dS > 

p p t^7r^(2p + l) 2 ' p \ L \ r i J <W> 2<W> 

Eq.(j21J) describes the contribution to the internal friction due to dislocation dipoles under 
the action of both periodic and random forces. As is known, the decrement describes how 
quickly the amplitude of a wave decays. In turn, the scattering rate, which is necessary 
to describe the thermal conductivity, defines the rate at which the incident wave loses its 
energy to the dislocation dipoles. It is given by 

= ^. (25) 

7T 

Notice that both the decrement and the scattering rate depend on two sets of parameters. 
The first one (T , B, m, D) characterizes the string and can be calculated from Eqs. (0) 
and (jllj) . The second set comes from the characteristics of the dipole's ensemble such as the 
density of dipoles A and the correlation function < 7](T — S)r](T) >. In general, the disloca- 
tion density depends on the plastic deformation in a sample. Typically, a phenomenological 
relation between the plastic strain e and the dislocation density is used. Experimentally 
the plastically deformed samples of Al, Ta, and Nb were studied . li2i ? i4i fr The densities of 
dislocations in these materials depend on the specimen preparation and vary in the range 
10 13 — 4 x 10 14 m" 2 . 1 ! 2 ! 3 ! 4 ! 5 In our calculations we try to choose such dipole densities which 
give the best fit. 

As was mentioned in the introduction, we assume that the stochastic stress component 
appears naturally in a high-density dislocation arrangement. Let us discuss this point in 
more detail. Evidently, it is impossible to describe the exact dynamical behavior of large 
dislocation ensemble in the complex potential relief. Instead, one can use the statistical con- 
sideration by introducing the statistical fluctuations of stresses. Notice that even idealized 



numerical simulations of the overdamped dynamics of an array of parallel edge dislocations 



with long-range interactions in Ref. |22j show an appearance in this system of the stochastic 
component of the internal stress field. According to Ref. |22j, for a period of time shorter than 
the relaxation time of the dislocation system, the internal stresses caused by dislocations in 
the arrangement can be viewed as a sum of slowly varying mean stress originating from the 
internal stress of the smoothed out dislocation distribution and a rapidly varying, highly 
irregular function with a zero mean value that satisfies the requirement of no correlation at 
different times or places and that represents the influence of the nearest neighbors. In this 
case, the stochastic stress component has a white noise character. 

Our calculations show that the white noise does not give a contribution to the internal 
friction. As is well known, however, the white noise concept does not reflect the real physical 
processes in complex systems. A possible way to generalize the consideration is to introduce 
the colored noise which can be modeled by exponentially decreasing correlator in the simplest 
case. As possible sources of colored noise in dislocation ensemble one can mention an external 
force, a long-range interaction between dislocations, the lattice itself, etc. Actually, we 
consider the dichotomic noise which is characterized by 

< rj(t) >= 0, < vitMh) >= vle- Mtl - t2 \ (26) 

where a is a parameter of the corresponding Poisson process, rj^ the dispersion. Taking 
into account Eqs.(j2U) an d (ESI), as well as the fact that for dichotomic noise < W >= 
(t}q + <To/2)/2/z, the total decrement reads 



with 



(Aa 2 m + 2aB + D) 



2\ Un/Lf 

Figffl shows the decrement given by Eq. (j27|) as a function of the reduced frequency for 
parameters of polycrystalline Al. The corresponding material constants for Al, Ta and Nb 
are presented in Table |l] As is seen from Fig^ at low frequencies (up to 10 8 Hz) the total 
decrement is a constant. It is important to note that such behavior is entirely determined 
by the contribution from the random part of stresses. The situation changes drastically 
near the resonance frequency where the decrement has a sharp peak and mainly depends 
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on losses due to the periodic external stresses. Namely, A p is found to be 10 4 larger than 
A r near the resonance frequency. Notice that a similar behavior is also found for both Ta 
and Nb. In experiments,—^ the low-frequency (OOkHa 1 ^ and 1kHz-) internal friction was 
studied. Therefore, the region of the constant decrement is of our main interest here. In 
addition, in our model the following condition is fulfilled 

D > T tt 2 /L 2 . (28) 

This is a variant of the so-called "rigid" relief approximation.— Moreover, a reasonable 
restriction on the dispersion of colored noise can be imposed in the form ^q/cJq -C 1. In view 
of Eq. (j28|) . performing the summation in Eq. (}2"Tj) one finally obtains 

87r/^ 2 A ?7 2 a 2nfiABb 2 u 
— p; 2 1 n2 • 

The decrement given by Eq. (|2T?|) for different specimens is presented in Table |n] As is seen, 
there is a good agreement with the experimental dat a 1 ! 4 ' 5 at appropriate choice of the model 
parameters. Notice that only two arbitrary model parameters (a/u> and ^o/^o) were used 
in our calculations. 

It was observed in experiments^ 4 * 5 - that the low-frequency internal friction does not de- 
pend on either the temperature or the frequency at low temperatures. In the framework 
of our model this behavior can be explained only when both the interaction between dis- 
locations (in a dipole) and the colored noise are present. Namely, the colored noise gives 
a dominant contribution (in comparison with the Granato-Lucke mechanism) to the decre- 
ment at low frequencies. At the same time, the athermal long-distance interaction (de- 
termined completely by the structure of the dislocation ensemble) provides the athermal 
behavior of the low-frequency internal friction at low temperatures. Indeed, this interac- 
tion ensures the validity of the "rigid" relief approximation (see Eq.(28)) when both the 
temperature-dependent terms (containing the damping parameter B) and the terms con- 
taining the tension and the mass of the string are negligible in Eq.(27). Let us note that 
a similar role of the long-range interactions was earlier discussed in Ref. |2J in a study of 
solid solutions. In particular, it was shown that the long-range interaction of dislocations 
with solute atoms distributed in the bulk of the crystal supports the athermal behavior of 
the nonlinear dislocation strain-amplitude-dependent internal friction in solid solutions . 24 ' 25 

We also assume the linear response of the dislocation array to the external periodic forces 
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which means a/c<j=const, r^/cr^ const. Without loss of generality one can fix the first con- 
stant as unity at numerical evaluations. Notice that, strongly speaking, these assumptions 
do not provide the independence of the decrement from frequency. Actually, this is true 
only for low-frequency asymptotic of the contribution to the decrement from random forces. 
Indeed, as is seen from Fig.l the decrement A r decreases rapidly at high frequencies. 



IV. THERMAL CONDUCTIVITY 



As is well-known, at low temperatures the main contributions to the thermal conductiv- 
ity come from the lattice defects and boundary scattering. Within the framework of the 
relaxation time description and the Debye approximation the lattice thermal conductivity 
is given by 




where the sum includes the longitudinal and two transverse acoustical phonon modes, u)Dj 
is the Debye frequency, Vj the sound velocity, Tj the total phonon relaxation time and Cj the 
specific heat of the phonon mode j. Assuming that all phonon modes are scattered almost 
equally by the dislocation dipoles, the total relaxation rate r~ l is written as 

r -i = r -i +r -i, (31) 

where r c _1 is the relaxation rate of the boundary scattering, and t^ 1 is the relaxation rate 
due to the dynamic phonon-dipole interaction. As another approximation, we will use the 
averaged sound velocity v instead of sound velocity of the phonon mode j. In other words, 
we drop the polarization effects. It is also convenient to use the dimensionless variable 
x = hu/kT so that Eq. ()3U|) takes the form 

e/T 

k = (k 3 /2Tr 2 vh 3 )T 3 J x 2 r(x)C(x)dx, (32) 
o 

where 6 is the Debye temperature. As is well-known, the boundary scattering gives a con- 
stant r^ 1 which is determined by geometry of samples. This contribution can be estimated 
by using the known Casimir formula, r^" 1 is the function of both the decrement A t and the 
frequency u in accordance with Eq.(|25j). The results of numerical calculations for all the 
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materials are shown in Figs. 2-4. The dislocation length is chosen to be 10~ 6 m in agreement 
with the acoustic experiments while the dipole separation I varies from one material to an- 
other to get the necessary resonance frequency. The best fits were found for the densities 
Aai = 6 x 10 12 m -2 and = 10 13 m~ 2 . The calculated curve for Ta is shown in Fig. 4 
with Ato, = 2 x 10 14 m~ 2 . Notice that these densities of dipoles differ from those used in 
the description of the internal friction experiments because the thermal conductivity and 
internal friction were actually measured on different samples.— i 2 i 3 i 4 i 5 Therefore, to clarify the 
problem it would be desirable to measure both the internal friction and the thermal con- 
ductivity on the same samples under the same conditions (to avoid stress relaxation and 
preserve the initial characteristics of dislocation ensembles). It is important to note that 
the influence of dichotomic noise on the thermal conductivity is found to be negligible. The 
reason is that the main contribution to the thermal conductivity comes from the region near 
the resonance frequency of optically vibrating dislocation dipoles. In this case, the term 
due to the colored noise can be disregarded in the scattering rate ( in Eq. (j25|) A r becomes 
much smaller than A p , see Fig.l). Thus, in contrast to the internal friction, the dichotomic 
noise does not influence the thermal conductivity. As is seen from FigsGKJ there is a good 
agreement with the experimental data at low temperatures. Notice that a rapid increase of 
experimental curves with temperature begins in the region where the electronic contribution 
becomes essential due to an increase of a number of normal electrons in superconducting 
samples. 



V. SUMMARY 



In this paper, we have shown that the presence of a stationary random process in disloca- 
tion ensemble caused by the external perturbation markedly modifies the decrement in the 
low-frequency range. First, the decrement does not depend on either frequency or temper- 
ature at low temperatures and, second, its value becomes much higher in comparison with 
the Granato-Lucke classical case of periodically affected dislocations. Both these findings 
agree well with the experimental data. At the same time, an explicit type and probabilistic 
characteristics of the proposed stochastic process are still not understood because dislocation 
ensembles show a very complicated dynamics. In fact, by introducing the dichotomic noise 
we simplify the problem of the nonlinear dynamics of dislocation dipoles to the linear one. 
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We have considered the influence of both external force and colored noise on the internal 
friction. We have also calculated the contribution to the thermal conductivity due to dislo- 
cation dipoles. The concept of dislocation dipoles allows us to obtain the required resonance 
frequency and, as a consequence, to describe the experimental results for all three plastically 
deformed metals. It was found that the dichotomic noise does not play any essential role in 
the thermal conductivity. 

Nevertheless, it is still impossible to say that the colored noise in the dislocation ensem- 
ble is the only way to explain the experimentally observed damping in Refs. QQ, and H 
The purpose of our paper was to draw attention to the fact that this mechanism gives an 
important contribution to the damping at low temperatures. A good agreement with the 
experiment was obtained by appropriate choice of two arbitrary parameters noticed above. 
There are other known mechanisms of damping which could be of importance to explain the 
experimental data. 

One of them was suggested in Ref. Q where the tunneling of some entities (presumably dis- 
locations or dislocation kinks) was discussed as a possible reason for a large internal friction 
observed in plastically deformed metals at low temperatures. Theoretically, the structures, 
energy barriers, effective masses, and quantum tunneling rates for dislocation kinks and jogs 
in copper screw dislocations were determined in Ref. |2y by using the molecular-dynamics 
calculations. It was found that for dislocation kinks in screw dislocations both the energy 
barrier and the effective mass are markedly reduced so that tunneling should occur readily. 
The kink tunneling is very sensitive to both the effective mass and the Peierls-like barrier 
for migration along the dislocation which, in turn, depend on the spatial extent of the kink. 
In particular, for wide kinks the barrier grows exponentially with decreasing kink width 26 . 
Hence, even a minor change of the kink width results in the substantial variation of the 
WKB factor {suppressing quantum tunneling). Notice that the value of the kink width used 
in Refs. l2ly m for Al, Ta and Nb (w = 5b ~ 1.4 nm) is about three times shorter than 
that in Ref. 2y for Cu. Besides, typically kink tunneling is suppressed by dissipation (see, 
e.g., Ref. l27l ). In addition, unfortunately both the resonance frequency and the scattering 
cross-section of the tunneling process are not yet available. This gives no way of calculating 
the contribution to the thermal conductivity within the kink tunneling model. Therefore, 
additional investigations within the tunneling kink model are necessary to describe a series 
of experiments in plastically deformed Al, Ta and Nb.-i* 2 ^^ 
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As another possible mechanism, the concept of kinks moving in the secondary Peierls 
potential was suggested for the description of observed experimental data.-^sM Despite a 
good agreement between the results of the kink model and the experiment, two comments 
are in order. First, it was found that a good fit for the thermal conductivity takes place 
only for a surprisingly large number of kinks (~ 500) per dislocation line of length 10 _6 m. 
In this case, the separation between kinks turns out to be 2nm. This value is comparable 
with the kink width w ~ 1.4 nm used in Refs. 00Q. At the same time, the radiation 
mechanism proposed in Ref. |28| is essentially based on the long-range interaction between 
kinks and is only valid when the separation between kinks is large compared to the kink 
width (otherwise the interaction between kinks may change drastically)^. Second, within 
the kink model the internal friction significantly (as a third power) depends on the Peierls 
stress whose reliable estimation is a complex experimental task^. As a final remark, at 
large densities the interaction between dislocations is of importance, hence the complex 
dislocation dynamics is to be expected. 
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TABLE I: Materials constant for Al (from Ref. 3), Ta (from Ref. 2), and Nb (from Ref. 4). 



Sample v fi (GPa) p (kg/m 3 ) c t (m/s) C/ (m/s) 



Al 


0.34 


24.8 


2690 


3290 


6600 


Ta 


0.35 


66 


16600 


2040 


4230 


Nb 


0.35 


47.3 


8560 


2140 


4720 
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TABLE II: Low frequency internal friction calculated according to Eq.(29) for different metals. 



Sample Plastic strain(%) A (m 2 ) 2na/uj r] /a a A t 

theory exp. 



Al» 


10 


10 14 


1 


1.75 x 


io- 2 


2xl0~ 4 


IO" 4 


Ta» 


15-20 


10 14 


1 


4.21 x 


io- 2 


1.5xl0~ 4 


1.5xl0" 4 


Nb c 


19 


10 14 


1 


9.935 x 


io- 2 


1.5xl0" 4 


1.5xHT 4 



"5N bulk Al from Ref. 1. 

b Ta 393c from Ref. 5. 

" Nb Z 48-250 from Ref. 4. 
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10~ 9 10" 7 10" 5 10" 3 10" 1 10 1 



FIG. 1: The calculated internal friction using the parameter set for the 10% deformed high purity 
5N polycrystalline aluminium sample (Tables I and II) and Qf = 10 11 Hz. 
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0.15 



FIG. 2: Thermal conductivity for Al. Experimental data denoted by squares are from Ref. 3. 
The solid line represents the theoretical curve for A^y = 6 x 10 12 vaT 2 using the parameter set from 
Table I. The resonance frequency for Al is flf = 10 11 Hz. 
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FIG. 3: Thermal conductivity for different Ta samples. Experimental data denoted by squares and 
circles are from Ref. 2. The solid line represents the theoretical curve for Ax a = 2 x 10 14 mT 2 
using the parameter set from Table I. The resonance frequency for Ta is £lf = 2 x 10 11 Hz. 
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FIG. 4: Thermal conductivity for Nb. Experimental data denoted by circles are from Ref. 4. 
Solid line represents the best fit with Am, = 10 13 mT 2 using the parameter set from Table I. The 
resonance frequency for Nb is flf = 5 x 10 11 Hz. 
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